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ABSTRACT 


A stability criterion is presented for spinning satel- 
lites, which consist of a rigid, symmetric, central body and 
four radial booms of arbitrary flexibility. 
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AniTUDE STABILITY OF A CLASS OF 
PARTIALLY FLEXIBLE SPINNING SATELLITES 


by 

Thomas W. Flatley 
Goddard Space Flight Center 

INTRODUCTION 

Spacecraft flexibility can invalidate the well known attitude stability criterion for spinning 
satellites. That is, rotation of a satellite about the maximum moment of inertia axis for the nom- 
inal spacecraft configuration may not be stable in the presence of energy dissipation. K struc- 
tural flexibility permits configuration changes and strain energy storage when the spacecraft spins 
about an axis other than the nominal spin axis, a new criterion is needed. 

The class of spacecraft considered here consists of a rigid, symmetric central body with four 
equally spaced uniform radial booms nominally in a plane through the center of mass and normal 
to the intended spin axis (Figure 1). 

STABILITY CRITERION 
Approach 

Consider a forced rotation of the system at 
a rate oj about a body fixed axis through the 
center of mass and at an angle a with respect to 
the nominal spin axis, with the booms in dy- 
namic equilibrium. Associate with this rotation 
a moment of inertia l(a, co) , a kinetic energy 
1/2 and a potential energy V(a, oS) , 

Consider also a rotation about the nominal spin axis with an equivalent angular momentum. 
Associate with this rotation an angular velocity a)^, a moment of inertia Iq, a kinetic energy 
1/2 Iq no potential energy. If the kinetic energy of this rotation is less than the total energy 

associated with the forced rotation, then rotation about the nominal axis is stable in the presence of 
energy dissipation. 

This stability criterion can be expressed as 

■2 + V , 



1 



2V 


or, 


which becomes 


2V 


for small a when i ~lo* 

Analysis of Forced Rotation 

Potential Energy 

Locate o> such that only one pair of booms is deflected. The potential energy of the system 
will be 


, Ji 

2EI 

•'o 


ds , 


where M = bending moment, Ei = boom stiffness factor, and I = boom length. 
K we let 

Ji - 1 - ^ 

EI " R " ds ’ 


the left hand side of the stability criterion becomes 


LHS 


A 

2EI I Y 

,2j Usj 


ds 


See Figure 2. 

Moment of Inertia 

The moment of inertia about the axis of forced 
rotation will be 


I - IjjCos^a + Ip sin^ a 


A 

+ -^/o[(a+^)^-a^] + 2 I pv‘^ ds , 

n 



Figure 2— Deflected boom geometry. 
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where Ijj and Ip are central body moments of inertia, a is the radius at which the boom is attached, 
and p is the boom's lineal mass density. 

Now, 


f-e 

^0 “ Ir + 3 p[(a + 2 1 p(a + s)^ ds 


SO 


rl 

I - Iq = (ip - Ip) sin^ a + 2p [r^ - (a + s)^] ds 


Now dr/ds = cos (a - 0) , and r(0) = a cos a. Thus for small a, 


r ^ all 






ds 


which becomes 


r (a + s) " " 2 “ 


a + 




ds 


Then 


(a + s)^ - a^(a+s) 


a + 


rti-)’ 


ds 


and 


2p 


I [r^-(a + s)^jds ~ -pa^’Sa(a+s)^ +1 2(a+s)j (— - l) 

•^0 0-^0 Jo 


ds 


ds 


Thus, for small a, the right hand side of the stability criterion becomes 


RHS = a* -^Ip - Ip -^at(2a +^) 


-,f L,,„ r (4-.r 

0 •'0 


ds 


ds - 
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Introducing the dimensionless quantities 


I-.. 




and 


2a 

e = X ’ 


the stability criterion becomes 


2EI r* 



1 r 


du > Ip - Ir - ^ 


( e + 2u) 1 du 

Jo 


and a stability boundary can be written 


AI = (^) ^ + + 




1 

1 

( € + 2u) du 

du 

-^0 

Jo 

J 


where 


pco 




2EI 


Series Representation 

Now, let 


It follows that 


N 

4> % Y' 4>„ u"-> 

n=l 




n=l j=l 


/*u N „ n 

l/’d- = 


^n+l-j ’ 


n=l j=l 


4 



I (€ + 2u) I cjp' du 
•'o Jo 


du 



n=l* j = l 


d0 


IM 


(n- u 


n~2 


r<-i 

E 

n = l 





N^l 


Lz: 

n-l j“l 




n-l 



N-l 


n=l j=l 


i) 


j+l ^n+2-j 


The stability boundary thus becomes 


AI - ^ ^ j(n + 1- ^ (1 + 6) 

L n=l j=l 

n=l j=l 

ROD ANALYSIS 

Equilibrium Beam Equations 

The equations of equilibrium for a thin rod and the constitutive relationship for an elastic rod 

are 



di^ ^ P - ^ ’ 


dM 

^+uxF + m - 0, 


and 


M = 



See Figure 3. 



♦Hildebrand, Francis B., "Advanced Calculus for Applications," Prentice-Hall, Inc., Englewood Cliffs, N. J., 1964. 
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For the planar problem at hand, 

F = Tu + Qn , where T = tension, Q = shear; 
p = centrifugal loading = pcJ^ r; 

M = Mb, where M = bending moment, b = uxn; 
m = applied couples = 0; 

1 

5 = , where e = slope of rod; 

du dn _ d0 

ds ~ ds ** ’ ds ds ** " 

Now, the magnitude of r is given by 

V 

r = a cos a + I cos(a-^)ds , 

Jo 

and a unit vector parallel to r is 

r = u cos (a - ^) + n sin (a - ^) . 

Thus keeping only first order terms in the small angles a and 0 we have the approximation 

r ~ (a + s) [u + (a - ^)n j . 

The equilibrium beam equations thus reduce to; 


dT di9 , 

di - Q (a + s) 

= 0 

d(9 

ds 

= 0 

dM ^ 
HI " Q 

= 0 

M 

d^ 

El 

ds 


The corresponding boimdary conditions are 



If we ignore the product of the small terms Q and 6.6/ As, integrate the first equation, and com- 
bine, we have 


e 

ds3 


_ 6 ^ 9 r „ ~ o T d0 

El - -J- [(a + '£.)2-(a+s)2] + pw2(a + s)(a-(9) , 

"" W {2a(a + s) + [2al + l^) ^ - 2a(s ^ +i5) - (s^ ^ +2si9)]- - 


0 ( 0 ) 


ds 


ds 


L/ . 

(^) 


0 . 


This equation integrates once, 


d^ e 

ds^ 



Since 


then 



0 , 


Now, making the substitutions 


c - “ a(a + . 


0 



s 


6 = 2^ 


we have 


d^ (/) _ ^2^4 v2 . .2') 

~ 2EI “'ll! +(^ + l)(<^+l) - eu(0+ 1) -U^ (4> + 1) - [2 + 1) J » 

d^ f 

- k{e + l-eu-u2}0 , 

= k{e(l -u) + (1 +u) (1 -u)} 0 , 

= k( 1 + e + u) ( 1 - u) . 
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The boundary conditions for this equation are 


4 ^( 0 ) = - 1 . 


d^(i) 


0 . 


Power Series Solution 

We now seek a power series solution of the form 


4> 



where N is large enough for convergence. 
To satisfy the boundary conditions, 


- -1 


and 



n= 1 




N- 1 

E 

n = 1 




0 . 


Now 




n= 1 


(n- 1) (n- 2) u' 


n ~ 3 


= 2<?b, + 


IN 

604 (n-l)(n-2)0„ 


= 203 + 604 u + 


N-4 


2 ^ (n + 3 )(n + 2 ) 0 „. 


n - 1 


and 


k(e + 1) 


N 

^ k( € + 1 ) ^ , 

n= 1 


= “k(€ + l) + k(€ + l)<^2^ 


N-2 



n - 1 
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“ keu<^ 


= -z: 

r^ = 1 




keu 


N~ 1 

z 


ke(j 6 


- ku^ <;6 


N 



n= 1 


Term by term comparison of the resulting equation yields 


and 


^3 


(n +3) (n+ 

‘^11 + 4 


k 

- -j (1 + e) ; 

[fi + (1 + €) (p^] ; 
k(l - k4>„ , 

(n + 3) (n + 2) [(1 ^^n + 2 ‘^nj ■ 


n=l,2, ■••,N-4; 


■^n = (n-l)(n-2) . 

n " 5, 6, • ■ • , N . 


Now each of these terms will be the sum of a constant and a linear term in i.e. 0. = a. 

+ b. Thus it is convenient to express each as a complex quantity (a b, ). Then, when multiply- 
ing terms by constants and forming linear combinations of terms, the constant and <t>^ dependent por- 
tions of the result are automatically obtained. 

With this convention, 


= (- 1 , 0 ) , 

^2 = ( 0 , 1 ) , 

4>^ = (-(1 + £)k/2, 0) , 

= (ek/6, (1+ Ok/6) . 


g 



is then determined from the boundary condition as follows: 


Y', = (A, B) ^ A + B.^2 = 0 > 



Values for all other terms are then calculated according to 

4>^ - real(0j + ^2 (^i)] • 

The previously mentioned "stability boundary" can be written as a function of k, e, and N: 


AI 


Y (1 + €) + 


N- 1 n 

n= 1 j = 1 




Limiting Cases 

In the extreme cases of infinite beam stiffness, where no deflection is possible, and vanishing 
beam stiffness, where the beam behaves as though hinged, no potential energy is associated with 
the forced rotation. The stability boundary is then 


AI 




r 


(e+2u) 



du . 


In the rigid rod case (k = 0), 


(9 = 0 , 0 = - 1 , 


SO 


AI € 

— ~ - ■^(l + e)+ I (6 + 2u)udu 


pi 


r 


“3 + € + 
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Thus for stability 


Al = Ip - Ir < e + ^] • 

Because the roll and pitch moments of inertia of the undistorted satellite are 

C = ^ ^ , 

A - 3" ^ , 

respectively, and because e = 2(a/^), then 

C = + 2pl^ ^ ^ ’ 

A = Ip + + 6 + ij . 

Ip - Ir < c - Ip - A + Ip . 

C > A , 

SO that the stability criterion reduces to the familiar ’’rigid body” criterion. 
In the hinged rod case (k = oo), 

0 = a y <^ = 0 


SO 


AI e 

= 2 ( 1-0 . 

Consider the possibility of a forced rota- 
tion, with hinged booms, being stable. Then, 
the torque on the central body caused by the 
centrifugal forces on the hinged booms (Figure 
4) must be just enough to maintain a constant 
angular velocity in body coordinates, i.e. 

CO = 0 , 

X ’ 


CO ^ CO sin a , 
y ’ 


u> 



co^ = CO cos a . 


Figure 4— Steady spin with hinged rods. 
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The torque on the central body is 


•'n 


T - - “ 2(a sin a) I p(a cos a + s) ds , 


” 2 apcj^ sin a a-^ cos a + * 


The Euler equations reduce to one, 



For small a, 

= P-e3 f [i + 2 f] , 

which is identical to the stability boundary above. 

RESULTS 

The ’^stability boundary’' dealt with above represents a maximum value of isi/pl^ for which 
the attitude motion of the spacecraft will be stable. K one plots this quantity as a function of the 
parameter k for some value of e , the region above the curve represents unstable combinations 
of variables and the region below stable configurations. Consideration of the limiting cases just 
described produces the values of asymptotically approached by the ends of the curve. For 

the general case o < k < 00 ^ points were calculated by the digital computation method described 
in this report. 

Figure 5 and Table 1 show the results obtained for intermediate values of k and several 
values of e. The dashed lines indicate asymptotic values for k = 0 and k = 00. 

Parameter definitions pertinent to use of the figure are: 

Al = maximum tolerable moment of inertia difference for the central body, 


AI 

AI 

pV 
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Table 1 


Results in Tabular Form. 


N . 6 

k N. 

0,0 

0.1 

0.2 

0.3 

0.4 

i 

0 

.6667 

.7717 

.8867 

1.0117 

1.1467 

0.1 

.6564 

.7599 

.8734 

,9967 

1,1300 

1 

.5787 

.6722 

.7751 

.8875 

1,0092 

10 

.3023 

.3719 

.4508 

.5393 

,6374 

100 

.0997 

.1594 

.2289 

,3083 

. 

.3974 

CO 

.0000 

.0550 

.1200 

,1950 

.2800 


From the monotonic nature of the curves, it is apparent that an increase in spin rate has a 
destabilizing effect, while an increase in boom stiffness (Ei) is stabilizing. The effect of ler^th 
and density changes are not so obvious, but, in the limiting cases of very small and very large k , 
it can be shown that increasing either one has a stabilizing effect. This is probably also true in 
general. 

In addition to the obvious use of these nondimensionalized results to check for attitude sta- 
bility (given some arbitrary set of design parameters), one could use them to generate design 
charts of various types for specific applications. For instance, one could display the stabilizing 
ability of a given type of boom (i-e., p and El specified), attached at a known radius a , by plotting 
Ai versus boom length for various spin rates. An envelope for such curves can be determined 
from the ’Timiting case" stability boundaries above, i.e. 



- p ^"3 + 23^^^ + 2 a ^^^ , 

Curves for various spin rates could be constructed by interpolation between data points 
generated as follows: 
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1. Select arbitrary points (1, Ai) within the envelope. 

2. Compute corresponding e and M/pl^ for each. 




I 


3. From Figure 5, determine values of k for each point. 

4. Calculate corresponding o)^s from the equation 


Figure 6 shows a qualitative example of 
such a chart. The upper and lower boundaries 
of the shaded region correspond to the cases 
a; = 0 and CO - CO ^ respectively. For a given Ai 
and spin rate, the minimum boom length re- 
quired for attitude stability can then be found 
from the plot as shown. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, September 9, 1968 
124-08-05-24-51 


2EIk 

■ 



Figure 6—Typical qualifatlve design chart. 


NASA-Langley, 1969 31 
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